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Isospin diffusion is probed as a function of the dissipated energy by studying two systems 
58j^;_l_58f^j and ^*Ni+l^^Au, over the incident energy range 52-74ylMeV. Experimental 
data are compared with the results of a microscopic transport model with two different 
parameterizations of the symmetry energy term. A better overall agreement between 
data and simulations is obtained when using a symmetry term with a potential part 
linearly increasing with nuclear density. The isospin equilibration time at 52 AVleY is 
estimated to 130±10 fm/c. 
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1. Introduction 

Cohisions between nuclei with different charge asymmetries may carry important 
information on the structure of the nuclear equation of state (EOS) symmetry term 
in density and temperature regions away from the normal value, that may be en- 
countered along the reaction path. For instance, the symmetry energy behaviour 
influences reaction processes, such as fragmentatio n, p re-equilibrium emission, N/Z 
equilibration between the two collisional partners. Among the sensitive observ- 
ables, in semi-periphera l col lisions, one can study isospin diffusion, either through 
isospin transport ratios ^^l^l gr by using the isotopic content of light particle and 
IMF emission and the as ymme try (N/Z) of the reconstructed quasi-projectiles (QP) 
and quasi-targets (QT).lSEl2l The degree of equilibration, that is related to the in- 
terplay between the reaction time and the typical time for isospin transport, can 
give information about important transport properties, such as drift and diffusion 
coefficients, and their relation with the density dependence of the symmetry energy. 

We have performed this kind of investigations by studying isospin transport 
effects on the reaction dynamics for two systems, with the same projectile, ^^Ni, 
and two different targets, ^^Ni and ^^^Au, at incident energies of 52AMeV and 
74A MeV.EEl The N/Z ratio of the two composite systems is N/Z=1.07 for Ni+Ni 
and N/Z=1.38 for Ni-f Au. This choice gives access to isospin effects in different 
conditions of charge (and mass) asymmetry and to their evolution with, the energy 
deposited into the system. In the symmetric Ni -t- Ni system isospin effects are 
essentially due to the pre-equilibrium emission. On the contrary, in the charge (and 
mass) asymmetric reactions, one can observe isospin transport between the two 
partners. The dependence of these mechanisms on the symmetry energy behaviour 
is discussed. 



2. Experiment 

2.1. Experimental details 

^^Ni (179 ^g/cm^) and ^^''Au (200 /ig/cm^) targets were bombarded with ^^Ni 
projectiles accelerated to 52 and 74AMeV by the GANIL facility. The charged 
products emitted in collisions were collected by the 47r detection array INDRA. ^ 
All elements were identified within one charge unit up to the projectile charge. 
Moreover H to Be isotopes were separated when their energy was high enough 
(above 3, 6, 8 MeV for p, d, t; 20-25 MeV for He isotopes; -60 MeV for Li and 
^^80 MeV for Be). In the following we shall call fragments the products for which 
only the atomic number is measured (Z>5). As the on-line trigger required four 
fired modules of the array, the off-line analysis only considered events in which four 
charged products were identified. ^ 
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2.2. Event selection 

After eliminating events in which less than 90% of the charge of the projectile was 
detected, we build subevents containing the ensemble of charged products which 
have a velocity higher than half the laboratory projectile velocity (i.e. forward emit- 
ted in the nucleon-nucleon (NN) frame). We call this ensemble "quasi-projectile" 
(QP) without prejudice on the equilibration of any degree of freedom of the system 
so defined. This choice eliminates the detection biases on the quasi-target (QT) side 
and avoids considering products coming from the QT. Finally to get homogeneous 
samples we ask that the total charge of the QP is comprised between 24 and 32. In 
all cases 1.3 to 2x10^ events are kept, amounting to 14-18% of the reaction cross 
sections. ^ 

2.3. The sorting variable 
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Fig. 1. Correlation between E^iss/Ec 
and two asy-EOS. From. El 



and the impact parameter for the four studied systems 



To follow the isospin diffusion as a function of the violence of the collision, we 
chose to sort the events as a function of the dissipated energy, calculated in a binary 
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hypothesis without mass transfer: 

Ediss = Ec.m. — -^Vrol^, (1) 

with fi the initial reduced mass, Ec.m. the available energy and 

Vrel=V^^{^X-^ (2) 
-^target 

Vqp is the QP velocity, reconstructed from the velocity of all the fragments it 
contains. Most of the QP events have only one fragment. For the Ni+Ni system, 
about 25-30% of the events have two or more fragments w hile fo r the Ni+Au system 
this percentage is smaller (~15%). It is demonstrated in 12113 that the velocity of 
the QP is a good parameter for following dissipation. The dissipated energy was 
calculated following the same procedure in the BNV simulations (see next section), 
and it is shown in fig [T] that there is a close correlation between Ediss/Ec.m. and 
the impact parameter. A complete description of the experimental findings can be 
found in. ^ 



3. The BNV transport model 

We follow the reaction dynamics solving the BNV transport equation, that describes 
the evolution of the one-body distribution function according to the nuclear mean- 
field and including the effects of two-body collisions. The test-particle prescription 
is adopted, using the TWINGO code with 50 test particles per nucleon. The 
main ingredients that enter this equation are the nuclear matter compressibility, 
the symmetry energy term and its density dependence and the nucleon-nucleon cross 
section. We take a soft isoscalar equation of state, with a compressibility modulus 
K =200 MeV, which is favoured e.g. from flow stud ies or from the confrontation 
data-dynamical simulations at intermediate energies. I ^'^ l ^'^ l 

Two different prescriptions for the behaviour of the symmetry energy are used 
to study the sensitivity of the results to the considered parameterization. Indeed 
one can write the symmetry energy part of the EOS as the sum of a kinetic term 
and a potential term (to be multiplied by [{N — Z)/A]'^), often approximated, for 
the sake of comparison, as: 

%^(p) = %^(^f/^ + %(^r (3) 

A 2 po 2 po 

where pa is the nuclear saturation density. The value of the 7 exponent, valid close 
to Po, determines whether the equation is "asy-stiff " (7 >1, potential term continu- 
ously increasing with p) or "asy-soft" (7 <1, potential term presenting a maximum 
between po and 2po). In this paper we chose as asy-stiff a potential symmetry 
term linearly increasing with nuclear density (7=1), and as asy-soft the SKM* 
param eter ization in which the potential symmetry term can be approximated with 
7=0.5. 
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The free nucleon-nucleon cross section with its angular, energy and isospin de- 
pendence was used. Indeed, as there is no consistency between the mean field and 
the residual interaction in transport codes, there is no reason to adopt an in-medium 
correction which may be "valid" with one specific mean field. For instance good 
reproductions of nuclear dyn amics above 40AMeV were obtained using the free 
nucleon-nucleon cross section. I ^^ l ^^ l 

For the two reactions, we have ran different impact parameters, from b=4 fm 
to b=10 fm for the Ni-f-Ni system and from b=4 fm to b=12 fm for the Ni-|-Au 
system. For each impact parameter 10 events were produced (one event represents 
already the mean trajectory of the reaction), for the two cases of symmetry energy 
parameterization. The detailed analysis of the results of the simulations can be 
found in. — 

4. Isospin diffusion 

The measure of isospin diffusion ideally requires the determination of the ratio 
N/Z for hot quasi-projectiles, just after QP and QT separate. Experimentally it 
is however not possible to fully reconstruct the primary QP because neutrons are 
not measured and the masses of the heavier fragments are not known. We chose 
to construct an isospin-dependent variable from the isotopically identified particles 
included in the QP: 

(< N > / < z >)cp = E E^--/ E E^-- (4) 

where Ni, and Py are respectively the numbers of neutrons and protons bound in 
particle v , v being d, t, ^He, ^He, <^He, ^Li, '^Li, ^Li, ^Li, ''Be, ^Be, ^°Be; free 
protons are excluded. N^vts is the number of events contained in the dissipated 
energy bin considered. The variable (N/Z)cp was calculated twice: first considering 
particles forward emitted in the nucleon-nucleon frame (Vparticio > Vjfj.Qj/2), and 
secondly keeping only particles forward emitted in the QP frame (Vparticie > Vgp). 
Indeed as the BNV simulation does not allow to identify mid-rapidity particles and 
light fragments, the theoretical value of (N/Z)cp is calculated only from the parti- 
cles evaporated by the QP, which we experimentally approximate as the particles 
forward emitted in the QP frame. 

4.1. Isospin diffusion in BNV simulations 

The evolution of the N/Z ratio of the quasi-projectiles as a function of Ediss/Ec.m. 
is reported in fig. [2] (N/Z)qp increases with the centrality of the collision for the 
two systems and the two beam energies. For the Ni-|-Ni system the variation of 
(N/Z)qp with centrality is small, and attributed to pre-equilibrium emission. Little 
dependence on the EOS appears at 52 ^MeV, while (N/Z)qp grows slightly higher 
at 74 AMeV for the asy-stiff case. Indeed, for this system with a small neutron ex- 
cess, more protons are emitted during the pre-equilibrium stage due to the coupled 
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Fig. 2. Upper curves in each panel: isospin ratio of hot QP, (N/Z)qp, vs Edigg/Ecm. obtained 
in the simulation with two asy-EOS (blue stars for asy-stiff and pink squares for asy-soft). The 
evolution of the variable (N/Z)cp obtained after de-exciting the QP with SIMON is represented 
by the lower curves (red stars and green squares). The lines correspond to linear fits. Adapted 
from. El 



effect of Coulomb repulsion and a less attractive symmetry potential for protons. 
This effect increases with the incident energy. On the contrary, the asy-soft case 
tends to emit more pre-equilibrium neutrons leading to a lower N/Z ratio. 

The evolution with centrality is much more pronounced for the neutron-rich and 
asymmetric Ni-|-Au system. In addition to pre-equilibrium effects, isospin transport 
takes place between the two partners of the collision, and increases with the violence 
of the collision. (N/Z)qp is mostly higher in the asy-soft than in the asy-stiff case 
for the two energies. The largest value reached, corresponding to b=4 fm, is lower 
at 74 ylMeV than at 52 ^MeV; this may be attributed to the shorter reaction 
times, and to the fact that the collision becomes more transparent. Thus the N/Z 
diffusion appears related to the degree of dissipation reached in the system and to 
the driving force provided by the symmetry term of the nuclear EOS, that speeds 
up the isospin equilibration among the reaction partners. An asy-soft EOS, more 
dissipative, favours isospin equilibration between the two partners, as found also in 
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other recent theoretical investigations. I " l *^ * l ^ * ^ ! 

4.2. Effects of secondary decay 

In figure [2] are also plotted the results concerning the variable (N/Z)cp, calculated 
after de-exciting the hot primary QPs with the help of the SIMON code. ■ The 
values of (N/Z)cp are always smaller and the evolution with dissipation is generally 
flatter than that of the N/Z of the primary QP: secondary decay weakens the isospin 
effects. 

At 52 ^MeV however the differences between the results of the two param- 
eterizations are more pronounced for (N/Z)cp, with respect to (N/Z)qp. Indeed 
excitation energies are larger in the asy-soft case, which favours the emission of 
neutron-richer particles, and thus enhances the effect due to the larger N/Z value 
observed in that case for the QP. At 74 A MeV, no discrimination between the asy- 
EOS is possible for Ni-|-Au (as was already the case for hot QPs), whereas for the 
reaction Ni-|-Ni the effects due to the de-excitation modify the initial trend imposed 
by the dynamical evolution: the values of (N/Z)qp are larger in the asy-stiff case; 
conversely the associated (N/Z)cp are smaller than those obtained in the asy-soft 
case. 

Finally (N/Z)cp, which appears in the simulations to be linearly correlated with 
the N/Z of the primary QP, is a good indicator of isospin transport effects and is 
sensitive to the asy-EOS. 

4.3. Comparison between experimental and simulated data 

The evolution of the experimental values of (N/Z)cp with dissipation is displayed 
in fig [3l Open points show the values obtained forward in the NN frame. Let us 
remind that in this case we mix mid-rapidity particles and those coming from the 
QP de-excitation. For the Ni-|-Ni system at both incident energies, (N/Z)cp varies 
by at most 1.5% when dissipation increases. This is the expected behaviour for 
this symmetric system where N/Z is only modified by pre-equilibrium emission, as 
explained in the previous subsection. For the Ni-|-Au system the isospin ratio is 
higher than that of the Ni-f Ni system whatever the dissipated energy. For the three 
first bins of dissipated energy (N/Z)cp has the same value for the two incident 
energies and slightly decreases with dissipation; this is due to the on-line trigger. 1^1 
The isospin ratio of the Ni-|-Au system is however higher than that of the Ni-|-Ni 
system, which could arise from the neutron skin of the Au target and/or from 
the mid-rapi dity pa. rticles included in our quasi-projectile selection which are more 
neutron rich. I ^^ l ^^ l 'pjjjg result is a first indication of isospin diffusion. At higher 
dissipated energies, (N/Z)cp presents a significant increase with dissipation and 
reaches higher values at 52^ MeV, while the trend is flatter at 74^ MeV. This 
may be interpreted as a progressive isospin diffusion when collisions become more 
central, in connection with a larger overlap of the reaction partners and thus a longer 
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Fig. 3. Isospin ratio of complex particles, (N/Z)cp, vs E(jiss/Ec.m. ■ Circles correspond to experi- 
mental data: open, forward in NN frame, close, forward in QP frame. Error bars are within the 
size of the symbols. Dotted and solid lines are as in fig [2] and hatched zones give error bars from 
simulations. Adapted from. El 



interaction time. For a given centrality, the separation time is longer at 52AMeV 
than at 74^ MeV, leaving more time to the two main partners to go towards isospin 
equihbration. 

The close points in fig [3] are related to the values of (N/Z)cp forward in the 
QP frame. They are in all cases smaller than the previous ones, and for Ni+Au at 
both energies, they grow faster with dissipation. This is because the mid-rapidity 
particles are no longer included: it is known that these particles are more neutron- 
rich, and that their isospin content is independent of the violence of the collision. 1^ 
The values of (N/Z)cp forward in the QP frame are compared with the results of 
the simulation, displayed in fig [3] by the lines and the hatched zones. A first result 
worth mentioning is that the chemical composition (N/Z) of the quasi-projectile 
forward emission appears as a very good representation of the composition of the 
entire quasi-projectile source. Such an observation seems to validate a posteriori 
the selection frequently used to characterize the QP de-excitation properties. 

When looking globally at the results for the four cases treated here, the agree- 
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ment is better when the asy-stifF EOS is used, i.e. a linear increase of the potential 
term of the symmetry energy around normal density. Note however that for Ni+Au 
at 52 ^MeV, where isospin transport effects are dominant, the close points lie in 
between the simulated results with the two EOS. This observation allows us to put 
an error bar on our result, expressed as 7=1±0.2 (see Eq. (jSj). This value is in 
reasonable agreement with that recently obtained from isospin diffusion in Sn+Sn 
systems at 50AMevS!, and matches the one derived from the competition between 
dissipative mechanisms for Ca+Ca,Ti at 25 A MeV. 

5. Isospin equilibration time 
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Fig. 4. Isospin ratio of primary Ni QP vs the 
impact parameter obtained in BNV simula- 
tions for the Ni+Au system at 52 A MeV. 
Adapted from. El 



The knowledge of the different time scales associated to the various degrees of 
freedom involved in heavy-ion collisions at intermediate energy is of crucial impor- 
tance to determine the physical properties of nuclear sources produced in the exit 
channel. Particularly with the announced exotic beams the N/Z degree of freedom 
will hopefully be explored over a wide range, and thus an estimate of the chemical 
(isospin) equilibration time becomes essential. It turns out that we can extract this 
information from the present data. Fig|4]shows, for the Ni-t-Au system at 52 A MeV, 
the evolution of the N/Z of the primary QP with the impact parameter obtained in 
BNV simulations. The equilibrium value, N/Z=1.38, is almost reached for b=4 fm, 
at the time when the QP and QT separate. This time is equal to 130±10 fm/c. A 
similar value of 115 fm/c has been reported for the Zn-|-Au system at 47AMeV. I^Sl 
Interestingly enough the equilibration of isospin can be derived directly from the 
experimental data. In the top-left panel of fig |3] one observes that the open and 
close points superimpose at high dissipation: the values of (N/Z)cp are the same 
at mid-rapidity and at velocities close to that of the QP. This is a strong indication 
of isospin equilibration. Note that more than 75% of the energy must have been 
dissipated before equilibrium is reached. Qualitatively this high value of the dis- 
sipated energy at isospin equilibration pleads for an asy-stifF EOS, as was shown 
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when con side ring; isospin transport ratios in simulated Sn+Sn collisions at 35 and 
SOAMeV.l^ni 

6. Conclusions 

We studied isospin transport as a function of dissipation or centrality in heavy-ion 
collisions for two beam energies, looking directly at the average isotopic content of 
the light particles emitted by the quasi-projectile. 

Experimentally we followed an isospin-dependent variable, (N/Z)cp, con- 
structed with the identified isotopes belonging to quasi-projectiles, as a function 
of the dissipated energy. Its value very slightly increases with centrality for the 
Ni-l-Ni system, for which no isospin transport effect is expected, whereas it evolves 
much more for Ni-|-Au, indicating the presence of isospin transport. 

Simulations show that for the Ni -|- Ni system, the N/Z of the quasi-projectile 
is essentially determined by proton rich pre-equilibrium emission, and thus slightly 
increases with centrality. The effect is more pronounced using an asy-stiff equation 
of state. For the Ni-|-Au system isospin transport takes place and the N/Z is larger in 
the asy-soft case. After de-excitation of the hot QPs, the isospin effects, although 
weaker, are still present and the variable (N/Z)cp is sensitive to the asy-EOS. 
The non-negligible effect of evaporation necessitates a precise knowledge of the 
evaporation properties of medium-mass exotic nuclei; the INDRA Collaboration 
has initiated an experimental program in that aim. 

The results of the simulations are analysed in the same way as the experimen- 
tal data. This demonstrates that (N/Z)cp is linearly related to (N/Z)qp. When 
considering the four cases under study we find a better overall agreement between 
simulations and experimental data for the asy-stiff case corresponding to a poten- 
tial symmetry term linearly increasing with nuclear density. Referring to Eq. ^ we 
propose a value 7=1 ±0.2. Information concerning the isospin equilibration time, is 
also obtained. At 52 AMeV for the Ni+Au and the most dissipative collisions we 
can infer from the data-model comparison that isospin equilibration is reached at 
130 ± 10 fm/c. 
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